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Nomenclature 


a,AOA 

angle-of-attack 

O^stall 

angle-of-attack at stall 

5 e 

elevator deflection 

AOA,a 

angle-of-attack 

c 

chord 

Cd 

section drag coefficient 

Cdo 

section drag coefficient at a=0° 

C h 

section hinge moment coefficient 

C l 

section lift coefficient 

Clitox 

peak section lift coefficient 

Cm,CMI/4 

section moment about 1/4 chord 

Cp 

pressure coefficient 

Ptotal 

total pressure 

Pstatic 

static pressure 

q,Q 

dynamic pressure 

V 

velocity 
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I. Summary 


In an effort to investigate the Ice Contaminated Tailplane Stall (ICTS) phenomenon, a full 
scale 2-Dimensional tailplane airfoil model of the NASA Glenn Research Center DeHavilland 
DHC-6 Twin Otter Icing Research Aircraft was tested in The Ohio State University 7x10 ft. wind 
tunnel. This test was a continuation of a similar wind tunnel test which took place in 1994. Two 
additional leading edge ice contours were added to the previous investigations’ test matrix. 

The primary Twin Otter tailplane consists of two airfoil sections. The first airfoil, 
Section 1, is used on approximately 70% of the entire tailplane span. This airfoil was tested 
uniced, and with four representative ice shapes attached to the leading edge, with angles of attack 
ranging from small positive to large negative values. The ice shapes were designated according to 
either the program that created them or their place of origin. The four ice representations were 
named S&C, LEWICE, IRT1 and IRT2. These shapes are representative of four potential ice 
formations. Aerodynamic coefficients, Cl, Cm, and Ch were obtained by integrating static 
pressure coefficient, Cp, values obtained from surface taps. Testing was accomplished for the 
most part at 60 kts, but some higher speed tests were run at 100 kts. The second airfoil. Section 
2, is used on 30% of the tailplane span and was also tested uniced, with maximum positive and 
negative elevator deflections. In addition to testing the representative ice shapes, a carborundum 
grit was applied to the leading edge to determine what, if any, effects were produced by a light 
debris buildup. 

Surprisingly, the generic ice shape models (S&C and LEWICE) produced more 
detrimental effects than the ice shape models molded from actual ice accretions. Perhaps this 
strange effect was caused by the jagged edge of the ice models causing turbulent boundary layers 
whereas the larger and more uniform S&C/LEWICE models caused intense vortices resulting in 
detached boundary layers. The effects of icing that would have the greatest effect on aircraft 
stability and control were shown in the changes in lift and hinge moment. For the clean 
Section 1 airfoil with S e = 0° elevator deflection, Qmax = —1.23 and (Xstaii = -17.8°. With the IRT1 
ice shape attached, these values were reduced to Qmax = -0.73 and ccs ta n = -10.5° and with the 
IRT2 ice shape attached, these values were reduced to Qmax = -0.85 and a sta ii = -15.5°. Hinge 
moment data with ice shapes showed a definite change with the LEWICE ice shape producing 
the largest shift. 

Velocity variations did not result in any noticeable differences with relatively small 
elevator deflections. Differences did show up once the elevator deflection reached the -20° 
position. Varying the airfoil Section did indeed have a significant effect on the aerodynamic 
coefficients. In general, the G max was slightly more negative for Section 2 than Section 1 while 
the difference in (Xstaii varied from zero to several degrees. 

The results obtained from this wind tunnel test have given an effective set of data with 
which to begin an analysis of tailplane aerodynamics in icing conditions. The data will be 
incorporated into the DHC-6 Twin Otter database and used to quantify tailplane aerodynamics 
during simulated maneuvering flight. This understanding will be needed and fully utilized in the 
development of maneuvers to identify aircraft susceptible to tailplane stall in icing conditions. 
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II. Introduction 


This wind tunnel test is the second test performed by The Ohio State University as part of 
the NASA, FAA, and OSU Tailplane Icing Program. The purpose of the program is to quantify 
the effect of tailplane icing on aircraft stability and control to aid in the development of a flight 
test procedure that will identify aircraft susceptible to tailplane icing. As a first step in this 
program, wind tunnel testing of the 2-Dimensional aerodynamics of the tailplane of the N ASA 
Glenn Research Center Icing Research Aircraft, a DeHavilland DHC-6 Twin Otter (Figure 1), 
was undertaken. Pressure data from surface taps were acquired and integrated to obtain lift, drag, 
pitching moment, and hinge moment coefficients. It should be noted that during the test of 1994 
(See Reference 1) pressure belt taps were applied to the surface of the model and used to collect 
pressure data, but evaluation and comparison of the data showed that the surface taps and the belt 
tap data were similar enough that the belt taps could be used on the flight vehicle to obtain valid 
pressures. The pressure belts were not used in the present tests. 

Two additional leading edge ice shapes were evaluated for the present test, as well as 
selected repetition of tests from the first test series of ice shapes designated S&C and Lewice. 
The shape of the two new leading edge ice accretion models, designated IRT1 and IRT2 
(Figure 5), were generated using the NASA IRT wind tunnel in an effort to provide actual ice 
accretion shape data before the next full scale flight test. The two new leading edge ice shapes 
are different from the previously tested S&C and Lewice (Figure 5) shapes in that they are actual 
casts of ice accreted on an airfoil; therefore having a relatively rough surface. The IRT1 shape is 
a model of an ice accretion collected for 22 minutes after de-icing boot failure whereas IRT2 is a 
model of ice after 15 minutes of accretion time with 3 minute cyclings of the 
de-icing boot. 

For this additional testing, the experimental facility as well as the model used were the 
same as in the initial test of 1994. This also included the traversing wake pressure measurement 
system. However, in an effort to keep up with time constraints, the Scanivalves were not used 
and a new system of pressure data acquisition by Pressure Systems Incorporated (PSr) was 
utilized. Using the PSI equipment reduced the average run time from 30 minutes to 
approximately 10 minutes. As noted earlier, the pressure belts were present on the model during 
the tests, but not used for surface pressure acquisition. It should be noted that the previous tests 
indicated that the tunnel speed of 100 kts did not significantly alter the aerodynamic coefficients; 
therefore the majority of data was taken at 60 knots. 

At the start of the present test sequence, a few of the earlier test points were re-examined 
in order to check for repeatability. After successful comparison with the 1994 series, new test 
points were added to expand the first tests and the new ice shapes were installed and tested. 
Lastly, a carborundum grit simulating a light icing encounter was applied to the leading edge and 
data was taken for comparison to other configurations. 
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Figure 1 NASA Glenn Research Center Icing Research Aircraft. 
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III. Experimental Facility 


The Aeronautical and Astronautical Research Laboratoiy Low Speed Wind Tunnel Facility 
(OSU LSWT) is a closed-loop circuit with a 7x10 ft. high speed test section, and a 14x16 ft. low 
speed test section. A 2000 Hp electric motor provides dynamic pressures up to 65 psf in the 7x10 ft. 
section (Figure 2). A sting balance mount allows models to be tested at angle-of-attack values 
ranging from -15 to +45 degrees, while the rotating tunnel floor has a range of -90 to 
+270 degrees. 1 



Figure 2 OSU Low Speed Wind Tunnel Facility 


Tunnel speed is manually set based on average values obtained from two dedicated tunnel 
total pressure probes and two static pressure probes. A PC-class computer is used to control sting 
and mounting table angles, and wake probe sweep, as well as data acquisition tasks. Up to 14 data 
channels can be sampled. Signal conditioning and amplification are provided by individual rack 
mounted units. 1 

Testing of 2-dimensional airfoils is conducted by installing a 7 ft. span model to fill the 
height of the high speed test section. The model is attached to the floor mounting table, which is 
then rotated to provide the desired angle-of-attack. Surface taps on the model are used to acquire 
pressure distributions, which are integrated to obtain model forces and coefficients. To determine 
momentum drag, wake survey probes are swept across the test section downstream of the model. 1 
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IV. Model Description 


a) Airfoil Sections 

The two airfoil sections tested are shown in Figure 3. The Section 1 airfoil is used on 70% 
of the tailplane span and is referred to as a clean airfoil when no ice shapes are attached to it. The 
tailplane area used by each section is indicated in the two- view aircraft drawing (Figure 1). 
Drawings of the four ice shapes are shown in Figure 5. Some minor smoothing of the lofted airfoil 
shape was required; a 0.020 inch "bump" at the lower leading edge was removed, and a straight line 
contour from the elevator main spar to the trailing edge was assumed. 1 

The primary factors in the model design were to provide a smooth aerodynamic surface 
while allowing flexibility in model configuration and low fabrication costs. A highly stiff design 
was also desired, to maintain model shape under load. The chosen model design consisted of airfoil 
sections machined in Machineable Plastic, also known as REN material, supported by metal tubing 
spars and aluminum ribs. This allowed the utilization of CNC machining to fabricate the airfoil 
shape with high accuracy. 1 

A schematic of the model is shown in Figure 4. The 20 inch span Machineable Plastic 
sections were supported by two internal metal tubing spars through each of the stabilizer and 
elevator sections. These sections transmitted their loads to the ribs through pins. The ribs contained 
keyways for keys mounted on the stabilizer spar tubes to provide stiffness under model moment 
loadings. Three ribs and two endplates structurally connected the spar tubes of the stabilizer to the 
forward spar tube of the elevator. The ribs of the elevator overlapped those of the stabilizer, to form 
a hinge around the forward elevator spar tube. End ribs on the elevator provided attach points to the 
endplates, which were used for changing the elevator deflection angle. The endplates contained 
rods to attach the model to the tunnel yaw table floor supports, and to a ceiling pivot. The model 
was mounted to the yaw table such that the table center of rotation was at 28% model chord. 1 

Removable plug sections were used to create the two airfoil sections tested. The stabilizer 
section consisted of the airfoil for Section 1, with an aft extension plug used to create the Section 2 
stabilizer airfoil. The elevator section consisted of the Section 2 airfoil with leading edge plugs to 
create the Section 1 airfoil. Figures 6 through 10 show the model as installed in the tunnel. 1 

Tap locations were chosen to provide high resolution where pressure gradients were large, 
and low resolution elsewhere to minimize the total number of taps. Taps were labeled starting with 
#1 on the leading edge of the stabilizer and ending with #86 on the elevator trailing edge. Odd 
numbered taps were located on the lower, or suction, side of the airfoil as installed in the aircraft, 
with even numbers on the upper, or pressure, side, except at the boat-tail ends. The tap locations of 
the ice shapes covered and replaced some of those of the clean airfoil. The final tap locations for 
each airfoil section and ice shape configuration are listed in Appendix A, which also shows a 
graphic of their positions on the airfoils. 1 

Surface taps were installed in-house prior to the completion of the section assembly. Taps 
were located on a 12 inch airfoil section using a Bridgeport milling machine with digital distance 
readout. Tap locations were determined relative to a coordinate system specifically chosen to be 
readily indicated using known positions on the section. Once the locations were marked along the 
section moldline, 0.0625 inch diameter holes were drilled normal to the surface 0.5 inches below 
one edge. Mating holes were then drilled vertically into the milled out end of the section. Stainless 
steel tubing was then fitted to these mating holes and routed along the section end and through 
tubing passageway holes. The tubing was glued and sealed to the mating hole using plumbers paste. 
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With the tubing completed, the milled out end was filled and glued to a mating 8 inch section to 
complete the 20 inch section. Strip-a-tube tubing was then attached to the stainless steel tubing and 
routed through the passageways in each section to connect the surface taps to the Scanivalve 
ports. The tubing length was approximately 15 ft. 1 

The belt was installed by first milling a slot into the pressure side of the stabilizer and 
elevator. Double sided tape and epoxy were then used to secure the belts to the section. The free 
ends of the belts were then routed through the slot and into the tubing passageways. 1 

b) Ice Shapes 

The representative ice shapes (See Figure 5) were named S&C, LEWICE, IRT 1 , and IRT2. 
The S&C ice shape has been used in several stability and control programs on the Twin Otter. This 
ice shape was determined using a combination of in-flight photographs of icing on the Twin Otter’s 
tailplane, and the ADS-4 report “Engineering Summary of Airframe Icing Technical Data”, by 
General Dynamics/Convair. The LEWICE ice shape was predicted using the LEWICE 1.3 
computer program. The specific input conditions were: velocity =120 knots, liquid water content = 
0.5 g/m 3 , median volumetric diameter = 20pm, icing time = 45 minutes, a = 0°, 
total temperature = -4°C. The ice shape was determined in five equal time steps. Each step 
consisted of a 9 minute accretion time, a redefining of the iced airfoil geometry, and a recalculation 
of the airflow over the new geometry. The 45 minute accretion time was chosen to be conservative 
while aligning with the FAA requirement of demonstrating capability with 45 minute accretion 
time on an unprotected surface 1 . 

The IRT1 and IRT2 ice shapes are models of actual leading edge ice accretion. The 
following method of manufacture was provided by Thomas Ratvasky of the NASA Glenn Research 
Center. First, ice is allowed to build up on the leading edge of an airfoil model in the NASA Glenn 
Icing Research wind tunnel. The ice located on either side of the area of interest is cleaned off of the 
model. Then a mold box with the bottom panel cut out in the shape of the leading edge is fitted 
against the model, forming an enclosed area surrounding the ice of interest and the leading edge of 
the airfoil model. The box is then filled with a silicon-based molding material (RTV 3110) that 
cures at below freezing temperatures. The molding material fills the box and surrounds the ice 
shape filling in all of the intricate cavities. (For more details about the tunnel conditions during the 
creation of IRT 1 and IRT2, please see Appendix F.) 

Once the molding material has cured, the mold box is removed from the airfoil model and 
the mold is extracted from the box. The mold has the shape of the box with a cavity that has the 
impression of the ice and leading edge of the model. 

Producing the mold is only the first of two steps in creating a finished ice shape. After the 
ice shaped mold is created, a plug of the airfoils’ leading edge is made and inserted into the ice 
shaped mold leaving a small cavity between the leading edge and the mold. A polyurethane 
material is mixed and poured into the cavity where it fills all of the ridges and bumps left in the ice 
mold. After curing, the casting is removed from the mold and plug. Because the properties of this 
type of mold are driven by the environment in the Icing Research Tunnel, they do not always have 
the resilience to separate nicely from the casting. In fact, some of the original mold may even 
tear off with the casting. These tom off pieces are cleaned off of the casting, which serves as the 
master casting. 
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The master casting is then used to make another mold using a different mold material that 
allows for a cleaner separation between mold and casting. This secondary mold is then used in the 
same manner described above to mass produce other castings. For this particular wind tunnel test, 
NASA made enough castings to span the airfoil model inside the 7x10 Low Speed Wind Tunnel. 
For the full scale flight test, there will be enough segments produced to fit on the entire tailspan 
(~20ft) of the DHC-6 Twin Otter. 

Because the installed ice shapes covered up many of the leading edge pressure taps, a 
method was devised to bring the pressure taps up through the ice shapes. This involved drilling new 
pressure taps and in some cases, simply attaching new tubing and rerouting the pressure lines. For 
the pressure taps in the vicinity of the de-icing boot, large holes were drilled through the ice shapes 
above where the taps were located. Wires were inserted through the large holes and into the 
pressure tap holes. Finally, the large holes were filled with body putty and sanded accordingly; the 
wires were then removed (See Figure 10). 

Lastly, a carborundum grit was applied to the leading edge to simulate insect debris on the 
leading edge. A template with a repeating “random” pattern of holes was laid over 6 inch wide 
double sided tape. The carborundum grit, type #30, was sprinkled onto the template with some of it 
sinking through the drilled holes and sticking to the double sided tape. This tape was then applied to 
the leading edge starting at approximately the stagnation line and laid down towards the trailing 
edge going over both the top and the bottom of the airfoil. A photograph of the grit tape is shown in 
Figure 1 1. A schematic of the applied grit tape is shown in Figure 12. 


Reference 1 (Baseline) 



Reference 2 



Figure 3 Section 1 and Section 2 airfoil profiles. 
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Figure 4 DHC-6 wind tunnel model schematic. 
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Figure 6 S&C ice shape installed. 


Figure 7 LEWICE ice shape installed. 
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Figure 11 Grit tape with grit applied. Shown actual size. 


Wing Section 
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V. Data Acquisition 


In an effort to reduce testing time, the previously used Scanivalves werenot utilized 
during this sequence of wind tunnel runs. Instead, a Pressure Systems Incorporated (PSI) data 
scanning system was used. This particular setup included a 780B Data Acquisition and Control 
Unit, a 780B Pressure Calibration Unit, an 81-IFC scanning module interface, one 32-port 5-psid 
pressure scanning module and two 2.5-psid 32-port pressure scanning modules. With the three 
pressure modules, a maximum of 96 individual pressures could be measured. The use of the 
PSI™ gear made it possible to complete a test run of 16 points in about 10 minutes whereas the 
Scanivalves™ would have taken approximately 30 minutes. 

The wake measurement system consisted of pressure probes mounted to a 16 ft. length of 
aerodynamic tubing which spanned the test section horizontally. The tubing was supported by 
roller guides located at each wall exterior and one roller guide mounted on a vertical section of 
aerodynamic tubing located at approximately 1/3 of the of the tunnel width from the port wall. 
Three sets of total and static pressure probes were attached to the aerodynamic tubing. These 
probes were 4 inches apart with 32 inches between sets to allow some overlap during the 36 inch 
sweep across the test section. An individual probe was used for each total and static pressure to 
allow measurements to be made at a constant tunnel cross-section location. Both total and static 
wake probe pressures were measured with 2.5 psi differential transducers. 1 

Tunnel conditions were measured using dedicated 2.0 psi differential transducers for pitot 
pressures while 0.5 psi differential transducers were used for the static pressures. All differential 
transducers were vented to the atmosphere with atmospheric pressure determined from a 
sensitive altimeter. 1 

A 386 class PC controlled all data gathering and data reduction operations with only 
tunnel airspeed being manually controlled. Angle-of-attack positioning, operation of the 
Scanivalves™, operation of the wake probe sweep, and recording of all raw data performed by the 
run program. 1 
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VI. Test Procedure 


Appendix B contains a run-log for this test. Repeatability runs of the Baseline 
configuration of 60 kts were taken at 8 e = 14.2°, 0.0°, -10.0°, and -26.6°. In addition, the last data 
point of each polar sweep was taken at a = 0.0° as a repeatability check. 

The daily test procedure consisted of first adjusting all signal conditioners and amplifiers. 
Then the atmospheric pressure was determined from a sensitive altimeter and recorded on paper 
for future insertion into the data acquisition program. The run software was set for the current 
test configuration and a no-wind tare was taken. The run program was then started, which 
positioned the model to the first angle-of-attack of the sweep (usually +4°). The tunnel fan was 
then started and manually set to the desired velocity as displayed by the run program. Once this 
velocity was stabilized, the angle-of-attack sweep mode was engaged. The run program then 
positioned the model and took data for all desired angle-of-attack points. A typical angle-of- 
attack sweep of 10 points (angles-of-attack) took approximately 10 minutes to complete. If a 
wake probe was desired, it occurred after pressure tap data were taken at each angle-of-attack. 
This sweep added an extra 2 minutes to the time taken for each individual test point. Because of 
this, only a limited number of wake probe sweeps were accomplished. An attempt was made 
between each test point to set the tunnel velocity at 60 kts; reason being that at high angles-of- 
attack, the rotating model created more blockage slowing the airspeed. 

With the angle-of-attack sweep completed, any desired post-run note was entered into the 
run program, and the raw data were saved to the hard disk. The raw data were the final result 
generated by the run program. The raw data file included the output voltages of all transducers, 
and the tap location and calibration data. This method of data storage allowed a run to be 
modified post-test if errors in the setup or tap-file were discovered after testing was completed. 

Limited tuft flow visualization was also performed. Upper and a lower rows of tufts 
consisting of white string were attached to the model with a black electrical tape. 
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VII. Data Reduction 


Aerodynamic coefficient data were obtained by integrating the pressure coefficient 
distributions of the stabilizer and elevator separately. A typical pressure distribution is shown in 
Figure 13. This integration was accomplished by applying the averaged pressure of adjacent taps 
to the area between them. When summed, these values gave coefficients defined relative to the 
local element coordinate system. These local coefficients were then converted to a lift and drag 
contribution by the appropriate axis transformation. Moment was obtained similarly using 
moment arms of the distance from a reference center to the center of each integration area. 

Wake data was only collected for a limited number of test runs. In order to adjust the 
pressure drag coefficient to include some wake drag data, a weighted-average bias value was 
determined by comparing the wake and pressure drag coefficients at a = 0° and 8 e = 0° for each 
section and ice shape. These bias values were then applied to all of the pressure drag coefficients 
obtained with the same section and ice shape configuration. 

Once the aerodynamic coefficients were obtained, a solid wall correction was applied. 
The formulas used to calculate the corrected data are shown in Appendix C. An example of the 
typical magnitude of the corrections is given in Figure 14 which shows corrected and uncorrected 
data for an ice shape run. A buoyancy term was not included in the corrections which is used in 
Ref. 2 as a C D correction. Ref. 2 implies that a buoyancy correction is valid only for large 
elements of an aircraft configuration and is not significant for airfoils alone. 1 

In addition to the aerodynamic coefficients, the reduction program also calculated other 
important data. Other output files generated contained uncorrected, corrected, wake drag, and 
wall static pressure coefficients. Appendix D contains a description of the file name and data 
structure conventions of these output files. 


DHC-6 Tailplane OSU 7'xlO' Wind Tunnel Data 
Run #2205, V=60kta, Section #2, a=-12.21 0 , 5e=14.2% Uncorrected Data 



Figure 13 Example of a typical pressure distribution. 
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Moment Coefficient, Cm . Lift coefficient. Cl 



Cl and C 


Section 2, 1RT1 Ice Shape, V=60kts, fie=(r 


j_ 



-10 -8 -6 

AOA (degrees) 


Comparison of Uncorrected Cm and Corrected Data 

Section 2, IRT1 Ice Shape, V=60kts, de=0® 


-16 -14 -12 -10 


-6 -4 -2 




AOA (degrees) 


Figure 14 Example of data reduction correction magnitudes 
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VIII. Results and Discussion 


The purpose of this test was to quantify the aerodynamic effects of various ice shapes on 
the twin Otter Tailplane. This quantification was sought relative to the tailplane stall 
phenomenon. The two most important phenomena were therefore lift and hinge moment. These 
subjects are discussed in the following section. Appendix E contains a complete data presentation 
of all test results of C L , C D , C M , and C H conveniently located on a single page for each. 

All data are presented with the tailplane as installed on the aircraft using an aircraft 
reference system; therefore the airfoil appears to be upside down. The values of interest are those 
occurring at negative angle-of-attack when the tailplane is generating negative lift, or positive 
down-force. Other sign conventions worth noting are; angle-of-attack positive with stabilizer 
leading edge up, moment and hinge moment positive when nose up, and elevator deflection 
positive with the trailing edge down. Figure 15 presents a schematic describing the axis system 
and sign conventions. 


4 Pitching 

Moment 4 Lift 


I 



- Lift 


4 Elevator 
Deflection 

+ Hinge 
Moment 


{ 4 Downforce) 


Figure 15 Schematic of axis systems and sign conventions. 


a) Ice Shape Effects 

The lift characteristics with zero elevator deflection are shown in Figure 16. It is clearly 
shown that the ice shapes reduce the Q max from that of the clean airfoil Ci max . Figure 16 also 
shows how the stalled angle-of-attack, a^u, is significantly reduced by the ice shapes. 
Specifically, the negative Ci m a X was reduced from -1.4 for the clean airfoil to approximately -0.7 
for the LEWICE ice shape configuration. The o^aii was also reduced from -16° to approximately 
—8° for similar configurations. Application of leading edge grit (the inclusion of grit data 
governed the use of this section 2 plot versus a section 1 plot) only affects the a^aii by a few 
degrees (-1° to 2°) and reduces (makes more positive) the Qmax by about 0.1. The hinge moment 
characteristics are shown in Figure 17. The Lewice and S&C ice shapes tend to have a more 
drastic affect on the hinge moment than the actual ice molds. 
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As may be expected, various elevator deflections cause the as, al i to either increase or 
decrease. This phenomenon is shown in Figure 18. Note that negative elevator deflections cause 
the (Xstaii to become more negative whereas positive deflections cause the a^u to grow more 
negative; much like how ailerons would affect an inverted wing. The exact opposite occurs with 
the Cl™* term. With the IRT1 ice shape installed (Fig. 18) and with 5 e = -10°, Ci max = 
-0.95 and a^n = -8.5°, whereas at 5 e = +10°, C ]max = -0.5 and 0^1 = -12.5°. These trends 
continue, although different in magnitude, even with large elevator deflections and occur with 
both Sections (Figs 19, 20, and 21). 

A similar effect happens with the hinge moment coefficient, C H (Fig. 22). A negative 
elevator deflection tends to raise the curve while positive deflections tend to lower the curve. The 
effect is even more noticeable with large elevator deflections (Fig. 23). With small deflections 
such as -10° and +10°, the difference between the curves is approximately ACh = 0.1. At large 
deflections, though, the difference approaches AC H = 0.13. This trend is continued and magnified 
with the use of Section 2 (Figs. 24 and 25). 

b) Velocity Effects 

The characteristic velocity effects are shown in Figures 26 through 29. These plots reflect 
the velocity effect on the IRT1 and IRT2 ice shapes along with the various sections. Figure 26, 
which shows the lift coefficient vs. angle-of-attack, demonstrates that as the elevator deflection 
grows more positive, the Ci max becomes less negative and oc ta ii becomes more negative. At 60kts, 
with 5 e = -20° there is a notable increase (more positive) in negative C L from approximately 3° to 
-T at a lower speed. With the IRT2 ice shape installed, the data (Fig. 28) show that there is even 
smaller or nonexistent difference between the high and low speed values. 

c) Repeatability 

In order to verify that comparison between the data of 1994 and 1996 would be valid, 
several repeat test runs were accomplished. These were carried out for both Sections 1 and 2 with 
elevator deflections of 14.2°, 0.0°, and -20.0°. The lift coefficients (Figs. 30 and 32) repeated 
relatively well except for the large deflection of -20.0° where the 1996 data showed a similar 
Qmax, but a more negative otstaii- The hinge moment data (Figs. 31 and 33) also showed such a 
general trend. The 1996 data matched fairly well except for the large deflection of -20.0°. 

d) Section Effects 

The Section 2 airfoil was tested at a velocity of 60 kts with 5 e = 14.2°, 0.0°, and -20.0°. Lift 
and hinge moment coefficient plots are shown in Figures 34 through 43. A notable characteristic of 
the Section 2 airfoil data is the general decrease (more negative) C L . For a 8 e of 0.0° and a clean 
airfoil, the difference between the Cimax’ s is approximately 0.15. This nominal difference is also 
seen with most of the other ice shapes, except for IRT2 (Figs 34 through 38). With the hinge 
moment data, there is also a difference between Section 1 and Section 2. Once again, the biggest 
differences exist for the elevator deflection of -20° (Figs. 39 through 43). 

As the elevator deflection grew more negative, the hinge moment coefficient for 
Section 2 became more positive than the Section 1 hinge moment coefficient. For the case of the 
elevator deflection becoming more positive, the hinge moment coefficient for Section 2 would 
become more negative than that of Section 1 . This trend continued throughout the hinge moment 
coefficient plots (Figs. 39 through 43). 
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Lift Coefficient, Cl 


DHC-6 Tailplane OdU 7X1 0 Wind Tunnel Data 
Ice Shape Effects, V=60 kts, Section 2, 6e=0.0 



Figure 16 Ice shape effects, lift coefficients. 


DHC — 6 Tailplane OSU 7X1 0 Wind Tunnel Data 
Ice Shape Effects, V=60 kts. Section 2, 6e=0.0 



Angle— of — Attac k (deg) 


Figure 17 Ice shape effects, hinge moment coefficients. 
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DHC-6Tai!plane OSU 7X1 0 Wind Tunnel Data 
Ice Shape Effects, V=60 kts. Section 1 , Min Deflections 



Figure 18 Ice shape effects, lift coefficient at small elevator deflections. Section 1. 


DHC— 6 Tailplane OSU 7X1 0 Wind Tunnel Data 
Ice Shape Effects, V=60 kts. Section 1 , Large Deflections 



Figure 19 Ice shape effects, lift coefficient at large elevator deflections. Section 1. 
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Lift Coefficient, Cl B Lift Coefficient, Cl 


DHC — 6Tailplane OSU 7X1 0 Wind Tunnel Data 
Ice Shape Effects, V=60 kts, Section 2, Min Deflections 



20 Ice shape effects, lift coefficient at small elevator deflections. Section 2. 


DHC — 6Tailplane OSU 7X10 Wind Tunnel Data 
Ice Shape Effects, V=60 kts, Section 2, Large Deflections 



Figure 21 Ice shape effects, lift coefficient at large elevator deflections. 
Section 2. 
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DHC — 6 Taiiplane OSU 7X1 0 Wind Tunnel Data 

Ice Shape Effects, V=60 kts, Section 1 , Min Deflections 


U- 4 - 


" 1 J 1 

-hj 1 1 r t rT ^ 

7 I i 

r*-u* r*i 


ana 


Clean, oe — — 1 0.0 
!RT1 , <5e=-1 0.0 
IRT2, <5e= — 1 0.0 
Clean, 5e = 0.0 
iRTI , 5e— 0.0 
IRT2, <5e=0,0 
Clean, <5e~1 0.0 (1 Q94 
iRTI . <5e = 1 0.0 
IRT2. 6e = 1 0.0 


LI LI 



-20 -18 -16 -14 -12 -10 -8 -6 -4 

Angle— of — Attac k (deg) 

Figure 22 Ice shape effects, hinge moment coefficient at small elevator 
deflections. Section 1. 


DHC — 6 Taiiplane OSU 7X1 0 Wind Tunnel Data 

Ice Shape Effects, V=60 kts. Section 1 , Large Deflections 






Clean, <5e= — 20 
IRTI , 6e= — 20.C 
IRT2, 5e= — 20. C 
Clean, <5e=0.0 
IRTI , 6e=0.0 
IRT2, <5e = 0.0 
Clean, 6e= 1 4.2 
IRTI , 5e = 1 4.2 
IRT2. 6e = 1 4.2 


-20 -18 -16 -14 -12 -10 


Angle — of— Attac k (deg) 

Figure 23 Ice shape effects, hinge moment coefficient at large elevator 
deflections. Section 1. 
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Hinge Moment Coefficient, Ch 


DHC-6 Tailplane OSU 7X1 0 Wind Tunnel Data 
Ice Shape Effects, V=60 kts, Section 2, Min Deflections 



Angle— of — Attac k (deg) 


Figure 24 Ice shape effects, hinge moment coefficient at small elevator 
deflections. Section 2. 


DHC-6 Tailplane OSU 7X1 0 Wind Tunnel Data 
Ice Shape Effects, V=60 kts, Section 2, Large Deflections 



Figure 25 Ice shape effects, hinge moment coefficient at large elevator 
deflections. Section 2. 


NAS A/CR— 2000-20992 1A/OL2 


23 




Hinge Moment Coefficient, Ch 


DHC — 6 Tailplane OSU 7X1 0 Wind Tunnel Data 

Velocity Effects, I RT 1 Ice Shape, Section 1 , Large <5e’s 



-20 - 18 - 16 - 14 - 12 - 10-8 -6 -4 -2 O 2 4 6 

Angle— of — Attac k (deg) 

Figure 26 Velocity effects, lift coefficient at large elevator deflections. 
Section 1. IRT1. 

DHC — 6 Tailplane OSU 7X1 0 Wind Tunnel Data 

Velocity Effects, IRT1 Ice Shape, Section 1 , Large die’s 



-20 -18 -16 -14 -12 -10 -8 -6 -4 -2 0 2 4 6 

Angle— of —Attac k (deg) 

Figure 27 Velocity effects, hinge moment coefficient at large elevator 
deflections. Section 1. IRT1. 
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DHC-6 Toil plane OSU 7X1 0 Wind Tunnel Data 

Velocity Effects, IRT2 Ice Shape, Section 1 , Large (5e’s 



Angle— of — Attac k (deg) 

Figure 28 Velocity effects, lift coefficient at large elevator deflections. 
Section 1. IRT2. 


DHC-6 Tailplane OSU 7X1 0 Wind Tunnel Data 
Velocity Effects, IRT2 Ice Shape, Section 1 , Large 6e’s 



Figure 29 Velocity effects, hinge moment coefficient at large elevator 
deflections. Section 1. IRT2. 
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Hinge Moment Coefficient, Ch Lift Coefficient, Cl 


DHC-6 Tailplane OSU 7X1 0 Wind Tunnel Data 

Repeatability, Section 1 , Clean, 60 kts 



Figure 30 Repeatability check, lift coefficient w/ 1994 data vs. 1996 
data. Section 1. 


DHC — 6 Tailplane OSU 7X1 0 Wind Tunnel Data 
Repeatability, Section 1 , Clean, 60 kts 



Figure 31 Repeatability check, hinge moment coefficient w/ 1994 data vs. 
1996 data. Section 1. 
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Hinge Moment Coefficient, Ch Lift Coefficient, Cl 


DHC-6 Tailplane OSU 7X1 0 Wind Tunnel Data 
Repeatability, Section 2, Clean, 60 kts 



Figure 32 Repeatability check, lift coefficient w/ 1994 data vs. 1996 
data. Section 2. 

DHC-6 Tailplane OSU 7X1 0 Wind Tunnel Data 
Repeatability, Section 2, Clean, 60 kts 



Figure 33 Repeatability check, hinge moment coefficient w/ 1994 
data vs. 1996 data. Section 2. 
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Lift Coefficient, Cl Lift Coefficient, Cl 


DHC-6Tai!p!ane OSU 7X1 0 Wind Tunnel Data 

Section effects, Clean, 60 kts 



-20 -18 -16 -14 -12 -10 -8 -6 -4 -2 0 2 4 6 


Angle— of —Attac k (deg) 

Figure 34 Section effects on the lift coefficient; clean. 


DHC — 6 Tailplane OSU 7X1 0 Wind Tunnel Data 

Section effects, S&C, 60 kts 



Figure 35 Section effects on the lift coefficient; S&C ice shape. 
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Lift Coefficient, Cl Lift Coefficient, Cl 


DHC-6 Tailplane OSU 7X1 0 Wind Tunnel Data 
Section effects, Lewice, 60 kts 



Figure 36 Section effects on the lift coefficient; LEWICE ice shape. 


DHC-6 Tailplane OSU 7X1 0 Wind Tunnel Data 
Section effects, IRT1 , 60 kts 



Figure 37 Section effects on the lift coefficient; IRT1 ice shape. 
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Hinge Moment Coefficient, Ch Lift Coefficient, Cl 


DHC— 6 Tailplane OSU 7X1 0 Wind Tunnel Data 
Section effects, IRT2, 60 kts 



Figure 38 Section effects on the lift coefficient; IRT2 ice shape. 


DHC — 6 Tailplane OSU 7X1 0 Wind Tunnel Data 

Section effects. Clean, 60 kts 



Figure 39 Section effects on the hinge moment coefficient; clean. 
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0 -8 


Angle— of— Attac k (deg) 

Figure 41 Section effects on the hinge moment coefficient; LEWICE ice 
shape. 
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Hinge Moment Coefficient, Ch Hinge Moment Coefficient, Ch 


DHC-6Tailplane OSU 7X1 0 Wind Tunnel Data 
Section effects, IRT1 , 60 kts 



Figure 42 Section effects on the hinge moment coefficient; 
IRT1 ice shape. 


DHC — 6Tailplane OSU 7X1 0 Wind Tunnel Data 
Section effects, IRT1 , 60 kts 



Figure 43 Section effects on the hinge moment coefficient; IRT2 ice 
shape. 
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IX. Conclusions 


A wind tunnel test has been completed which quantified the 2-D aerodynamics of the 
DHC-6 Twin Otter Tailplane with and without representative ice shapes. The results obtained are 
consistent with the expected aerodynamic trends. The effects of icing that would have the 
greatest effect on aircraft stability and control were shown in the changes in lift and hinge 
moment. For the clean Section 1 airfoil with 5 e = 0° elevator deflection, Ci max = —1.23 and oCstaii = 
-17.8°. With the IRT1 ice shape attached, these values were reduced to Ci max = - 0.73 and (Xs ta ii = 
-10.5°. For the IRT2 ice shape attached, these values were reduced to Ci max = - 0.85 and OCstaii = 
-15.5°. Hinge moment data with ice shapes showed a definite change with the LEWICE ice 
shape producing the largest shift. 

Velocity variations did not result in any noticeable characteristic differences and 
relatively small elevator deflections. Differences did show up once the elevator deflection 
reached the -20° position. 

Varying the airfoil Section did indeed have a significant effect on the aerodynamic 
coefficients. In general, the Ci max was slightly more negative for Section 2 than Section 1 while 
the difference in ots ta ]i varied from zero to several degrees. 

This wind tunnel test has produced an effective set of data with which to begin an 
analysis of tailplane aerodynamics in icing conditions. The data will allow the tailplane 
aerodynamics during simulated maneuvering flight to be quantified. Such information will allow 
a more complete understanding of the contribution of the tailplane to the aircraft’s motion during 
maneuvering flight. This understanding will be needed and fully utilized in the development of 
maneuvers to identify aircraft susceptible to tailplane stall in icing conditions. 
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Appendix A 
Tap Locations 


Section 1, Clean 


TAP# X/C 


1 

-.0011 

2 

.0003 

3 

.0001 

4 

.0041 

5 

.0039 

6 

.0100 

7 

.0100 

8 

.0175 

9 

.0169 

10 

.0250 

11 

.0250 

12 

.0373 

13 

.0380 

14 

.0500 

15 

.0500 

16 

.0750 

17 

.0750 

18 

.1000 

19 

.1000 

20 

.1256 

21 

.1260 

22 

.1500 

23 

.1500 

24 

.1753 

25 

.1750 

26 

.2000 

27 

.2000 

28 

.2500 

29 

.2500 

30 

.3000 


TAP# X/C 

31 .3000 

32 .3500 

33 .3500 

34 .4000 

35 .4000 

36 .4500 

37 .4500 

38 .5000 

39 .5000 

40 .5226 

41 .5243 

42 .5180 

43 .5185 

52 .5261 

53 .5363 

54 .5354 

55 .5462 

56 .5449 

57 .5759 

58 .5741 

59 .5955 

60 .5936 

61 .6152 

62 .6131 

63 .6346 

64 .6326 

70 .6569 

71 .6589 

72 .7009 

73 .7026 


TAP # X/C 

74 .7448 

75 .7463 

76 .7888 

77 .7901 

78 .8327 

79 .8338 

80 .8767 

81 .8775 

82 .9207 

83 .9212 

84 .9646 

85 .9650 

86 1.0000 


Notes: 

1) Stabilizer taps are #1^43, elevator taps are #52-86 
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Section 2, Clean 


TAP# X/C 

1 -.0011 

2 .0003 

3 .0001 

4 .0041 

5 .0039 

6 .0100 

7 .0100 

8 .0175 

9 .0169 

10 .0250 

11 .0250 

12 .0373 

13 .0380 

14 .0500 

15 .0500 

16 .0750 

17 .0750 

18 .1000 

19 .1000 

20 .1256 

21 .1260 

22 .1500 

23 .1500 

24 .1753 

25 .1750 

26 .2000 

27 .2000 

28 .2500 

29 .2500 

30 .3000 


TAP# X/C 

31 .3000 

32 .3500 

33 .3500 

34 .4000 

35 .4000 

36 .4500 

37 .4500 

38 .5000 

39 .5000 

44 .5500 

45 .5500 

46 .6000 

47 .6000 

48 .6173 

49 .6086 

50 .6084 

51 .6089 

65 .6119 

66 .6164 

67 .6174 

68 .6297 

69 .6314 

70 .6569 

71 .6589 

72 .7009 

73 .7026 

74 .7448 

75 .7463 

76 .7888 

77 .7901 


TAP# X/C 

78 .8327 

79 .8338 

80 .8767 

81 .8775 

82 .9207 

83 .9212 

84 .9646 

85 .9650 

86 1.0000 


Notes: 

1) Stabilizer taps are #1-51, elevator taps are #65-86 
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S&C Ice Shape Tap Locations 


TAP# 

X/C 

Y/C 

1 

-0.0084 

0.0030 

2 

-0.0027 

0.0201 

3 

-0.0014 

-0.0150 

4 

-0.0018 

0.0315 

5 

-0.0019 

-0.0264 

6 

0.0088 

0.0308 

7 

0.0088 

-0.0275 

8 

0.0175 

0.0286 

9 

0 . 0175 

-0.0258 


IRT1 Ice Shape Tap Locations 

TAP# X/C Y/C 


1 

-0.0118 

0 . 0128 

2 

-0.0098 

0.0258 

3 

-0.0214 

-0 . 0047 

4 

0.0026 

0.0226 

5 

-0.0112 

-0.0114 

6 

0.0084 

0.0253 

7 

0.0077 

-0.0189 

8 

0.0173 

0.0278 

9 

0.0152 

-0.0228 

10 

0.0250 

0.0285 

11 

0.0237 

-0.0266 

12 

0.0377 

0.0311 

13 

0.0365 

-0.0314 

14 

0.0502 

0 . 0325 

15 

0.0486 

-0.0358 

16 

0.0756 

0.0343 

17 

0.0738 

-0.0434 


LEWICE Ice Shape Tap Locations 


TAP# 

X/C 

i 

Y/C 

1 

-0.0104 

0.0029 

2 

-0.0207 

0.0338 

3 

-0.0165 

-0.0282 

4 

-0.0044 

0.0239 

5 

-0.0044 

-0 . 0213 


IRT2 Ice Shape Tap Locations 


TAP# 

X/C 

Y/C 

1 

-0.0057 

0.0033 

2 

-0.0037 

0.0110 

3 

-0.0055 

-0.0060 

4 

0.0005 

0.0176 

5 

-0.0009 

-0.0133 

6 

0.0095 

0 . 0238 

7 

0.0075 

-0.0189 

8 

0.0172 

0.0275 

9 

0.0148 

-0.0220 

10 

0.0253 

0.0290 

11 

0.0234 

-0.0257 

12 

0.0380 

0.0311 

13 

0.0367 

-0.0308 

14 

0.0506 

0.0322 

15 

0.0487 

-0.0353 

16 

0.0759 

0.0336 

17 

0.0741 

-0.0422 


Tap Locations for the Various Ice Shapes 
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-2 -1 01234 -2 -1 01234 

X(inche8) X(lnches) 


Airfoil Profiles With Ice Shapes and Tap Locations 
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Appendix B 
Run Log 


Run # Section 



Ice Shape 


C 


Clean 


Elevator 


mi 


Wake? 


No 




.HIS file 


trv.his 


Notes 



\wmm 


1 Lewice 


MB 

ce 


Lewice 



60 kts 


60 kts 


60 kts 


0 kts 






tip202.his 


tip204.his 








Bad run 


Bad run 


Bad run 


2229 

1 

Clean 

-2 



Bad run 


21 5. his I Bad run 


tip220.his 
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0060 1 

2 

S&C 

1 14.2°l 

60 kts 

no 

tiD237.his 

Repeat of 0010 


2 

S&C 

■KfjSi 

60 kts 

no 

tip238.his 

Repeat of 001 1 

0062 

2 

Lewice 

MBItBI 

60 kts 

ves 

tiD235.his 

Repeat of 0012 

0063 

2 ; 

Lewice 

■RSI 

60 kts ■ 

no 

tio234.his 

Repeat of 001 3 

0064 

2 

Lewice 


60 kts 

no 

tip236.his 

Repeat of 0014 
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Appendix C 

Solid Wall Correction Calculations 


The effect of the walls on the 2-Dimensional airfoil data must be considered. Corrections 
have been applied to the raw wind tunnel data using the methods of Ref. 2. For this test, three 
phenomena have been taken into account: 

Solid Blockage 

Solid blockage creates an increase in velocity due to the reduction of the area through 
which the air must flow. To quantify this effect, computations have been made by considering 
the 2-Dimensional model as a cylinder. This cylinder can be simulated as a doublet of strength 
p = 27tVa 2 (a: the cylinder radius), with the walls represented by the streamlines created by 
matching doublets on each side of the cylinder. Therefore, 

AV_a z 

K ~ K 

(The subscript u indicates uncorrected data) 


After doing a summation, we obtain. 


with: 


E 


sb 


fAV^l 


Act 


h: tunnel height 
c: model chord 


a - 


vv 

48, 


1 



x, y: airfoil coordinates, P is no-camber, symmetrical, pressure 
distribution. 


Usually A is given by a graph. In our case, A = 0.22 
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Wake Blockage 

Wake blockage creates a velocity increment at the model because of a pressure gradient 
due to higher velocity which keeps the flow around the model. In computations, the wake is 
simulated by a line source and the walls by an infinite vertical row of source/sink combinations. 
This reduces to: 

_ AV __£^ r 

“* V u 4 h du 

We can then correct the tunnel conditions by these equations: 

V = V.(1 + C) 

with : £ = e sb +£ wb 


q = q u {\ + 2e) 

Re = Re 1( (l + e) 

From the dynamic pressure effect and the wake gradient term, we get: 

C, = Cj l-3e„-2ej 


Streamline Curvature 

The lift and moment about the quarter chord of an airfoil are too large at a given angle of 
attack (which is also too large). This is due to the fact that the airfoil seems to have more camber 
because of the floor and ceiling. 

Calculations of this effect are made by assuming that the airfoil can be approximated by a 
single vortex at its quarter-chord point. The floor and ceiling are represented by a vertical row of 
vortices, extending to infinity and with alternating signs. The load on the airfoil can be 
decomposed as a flat plate loading, computed as an angle of attack correction and an elliptical 
loading which gives the lift, pitching moment, and hinge moment corrections. 

It can be shown that the upwash induced at the half chord by the two images is: 


A a = 


8/r 


J 

c~ 



Cl 
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We can consider that (c/4) 2 is small compared to h 2 , so it gives a equal to <7 previously found, 
and we have then : 


« = «« 


57.3(7 
+ 

In 


( 


C„+4C 


V 


\ 

/ 


C, = c, u ( l-(7-2e) 

aC, 

C , =C , (1 - 2e) + —4 

m— m—u 4 

4 4 


Hinge Moment 

To establish a correction for the hinge moment, we follow the same analysis process 
discussed above, but we consider only the flap or elevator instead of the entire airfoil. We are 
now at a three-quarter-chord point, so we lose accuracy by taking (3c/h) «h' in the Aa 
expression. Therefore, including this 3c/4 term results in o’=0.9c. 

So, 

c t =c„ a-2e)+jC f 

with: Cif = lift coefficient of the elevator. 
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Appendix D 

Data Reduction Output File Formats 


Uncorrected Coefficients 

File Name: AA500XX.DAT 
Run: OOXX 


Column Value 


1 

Angle-of-Attack (degrees) 

2 

Elevator Deflection (degrees) 

3 

Mach Number 

4 

Reynolds Number (millions) 

5 

Q (psi) 

6 

Velocity (kts.) 


From Surface Taps 

7 

CL (total) 

8 

CD (pressure only) 

9 

CM ( 1 /4c) 

10 

CL (Elevator only) 

11 

CH (Elevator) 


From Belt Taps (although not used) 

12 

CL 

13 

CD (pressure only) 

14 

CM (l/4c) 

15 

CL (Elevator only) 

16 

CH (Elevator) 

17 

Run Number 
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Corrected Coefficients using Surface Taps 

File Name: SC500XX.DAT 
Run: OOXX 

Column Value 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


Angle-of-Attack (degrees) 
Elevator Deflection (degrees) 
Mach Number 
Reynolds Number (millions) 
Q (psi) 

Velocity (kts.) 

CL (total) 

CD (pressure only) 

CM (l/4c) 

CH (Elevator) 
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Appendix E 

DHC-6 Twin Otter Tailplane Coefficient Data 


List of Figures From Appendix E. 

Figure El. Ice Shape Effects, Section 1, 5e = 0.0° 50 

Figure E 2. Ice Shape Effects, Section 1, 5e = 10.0° 51 

Figure E 3. Ice Shape Effects, Section 1, de= 14.2°. .................. ...... 52 

Figure E 4. Ice Shape Effects, Section 1, 5e=-10.0° 53 

Figure E 5. Ice Shape Effects, Section 1, 5e=-20.0° 54 

Figure E 6. Ice Shape Effects, Section 1, 5e=-26.6°........................................... 55 
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Figure E 2. Ice Shape Effects, Section 1, 5e = 10.0' 




Figure E 3. Ice Shape Effects, Section 1, de=14.2°. 








ce Shape Effects, V=60 kts, Section 1 , 6e 



Figure E 4. Ice Shape Effects, Section 1, 5e=-10.0°. 
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Figure E 5. Ice Shape Effects, Section 1, 5e=-20.0°. 
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Figure E 6. Ice Shape Effects, Section 1, 8e=-26.6' 
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Figure E 7. Ice Shape Effects, Section 2, Se=0.0°. 
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Figure Ell. Ice Shape Effects, Section 2, 6e=-20.0°. 
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Figure E 12. Section Effects, Clean, Large Elevator Deflections. 
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Figure E 13. Section Effects, IRT1 Ice Shape, Small Elevator Deflections. 








Figure E 14. Section Effects, IRT1 Ice Shape, Large Elevator Deflections. 












Figure E 15. Section Effects, IRT2 Ice Shape, Small Elevator Deflections. 
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Figure E 16. Section Effects, IRT2 Ice Shape, Large Elevator Deflections. 





Section Effects, LEWICE Ice Shape, V=60 kts, Large he’s 



Figure E 17. Section Effects, LEWICE Ice Shape, Large Elevator Deflections. 










Section Effects, S&C Ice Shape, V=60 kts, Large 6e’s 



Figure E 18. Section Effects, S&C Ice Shape, Large Elevator Deflections. 










Velocity Effects, IRT1, Section 1 
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Figure E 19. Velocity Effects, IRT1 Ice Shape, Section 1, 
Large Elevator Deflections. 
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Figure E 20. Velocity Effects, IRT2 Ice Shape, Section 1, 
Large Elevator Deflections. 






Clean, V= 60 kts, Section 



Figure E 21. Clean, Section 1 










Figure E 22. Clean, Section 2. 












Grit, V=60 kts, Section 2 
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Figure E 23. Grit Applied, Section 2 




IRT1 Ic e Shape, V=60 kts, Sec tion 1 
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Figure E 24. IRT1 Ice Shape, Section 1, V=60 kts. 





RT1 Ice Shape, V=1 00 kts, Section 1 



Figure E 25. IRT1 Ice Shape, Section 1, V-100 kts. 







» 

1.0 
0.5 
0.0 
-0.5 
-1.0 
-1 .5 

O 1 
CN 
1 

n cn t- 

d d 6 

O <N * 

odd 
1 1 

* 

lO'juapyjaooHn 


W0 , pJ°90t'/lP1U3! 

DIJ.j.^OQ }USlUO^ Buiqo^d 


NASA/CR— 2000-20992 1/VOL2 


75 



Figure E 26. IRT1 Ice Shape, Section 2. 







IRT2 Ice Shape, V=60 kts, Section 1 



Figure E 27. IRT2 Ice Shape, Section 1, V=60 kts. 








RT2 Ice Shape, V = 1 00 kts, Section 



10 '^uapupooun 


wq 'pjoqo t / 1 P psiDijpoo }usuuo^ 6u|n o>|d 


NASA/CR— 2000-20992 1/VOL2 


77 


Figure E 28. IRT2 Ice Shape, Section 1, V=100 kts. 










RT2 Ice Shape, V=60 kts, Section 2 



Figure E 29. IRT2 Ice Shape, Section 2. 







Appendix F 

IRT Ice Shape Information 


The information shown below was provided by Thomas Rat\’asky, NASA Flight 
Research Engineer. 


IRT1 Ice Shape 

Represents the leading edge ice buildup after failure of a de-icing boot. 


Airspeed: 
Angle-of-attack: 
Icing time: 
LWC: 

Total Temp: 
MVD: 


135 kts 
-2.9° 

22 minutes 
0.5 g/m 3 
24.3 °F 
20 microns 


IRT2 Ice Shape 

Represents the amount of leading edge ice accreted between de-icing boot cycles. 


Airspeed: 
Angle-of-attack: 
Icing time: 


LWC: 

Total Temp: 
MVD: 


135 kts 
-2.9° 

15 minutes with de-icing boot cycles every 3 minutes, so the final 
accretion is the residual ice left after the deicer cycle at 1 2 minutes 
plus the last 3 minutes of growth. 

0.5 g/m 3 
24.3 °F 
20 microns 
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